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A BIHARMDNIC RELAXATION METHOD FOR CALCULATING THERMAL 
STRESS IN COOLED IRREGULAR CYLINDERS 
By Arthur G. Holms 


SUMMARY 

A numerical method was developed for calculating thermal stresses 
in irregular cyl inder s cooled by one or more internal passages . The 
use of relaxation, methods and elementary methods of finite differences 
■was fo und to give approximations to the correct values -when compared 
with previously known solutions for concentric circular cylinders 
possessing symmetrical and asymmetrical temperature distributions . 


INTRODUCTION 

Use of Cooled Irregular Cylinders 

The evolution of aircraft propulsion systems has led to the fre- 
quent employment of cooled structures. The concentric hollow cylinder 
is a familiar example, although in some cases irregular cylinders such 
as cooled turbine blades with several internal passages (reference l) * 
are under consideration. A method of calculating stresses in thin- 
walled turbine blades of the air-cooled type is presented in refer- 
ence 2, but the general problem of calculating thermal stresses in 
long, hollow, thick-walled irregular cylinders has not been solved. 


Previous Work on Thermal Stresses in Hollow Cylinders 

Several methods of calculating thermal stresses for various special 
distributions of temperature in long hollow cylinders of particular 
shape have been developed. Some theoretical aspects of the general 
problem have also been discussed. In reference 3 the problem is 
regarded as an or din ary stress problem with given body and surface 
forces replacing the effects due to temperature distribution, whereas 
in reference 4 the equilibrium and boundary conditions of the theory 
of elasticity are used without modifications to exhibit the temperature 
effects as body and surface forces . The method of reference 4 was 
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applied to several special problems that had already been solved as 
well as to problems of* composite bodies and eccentric circular 
cylinders. No application of* analytical or numerical methods of 
calculating thermal stresses in cooled irregular multiply connected 
cylinders has been published. 


Scope of Present Investigation 

An investigation was conducted at the NACA Lewis laboratory to 
calculate thermal stresses in cooled irregular multiply connected 
cylinders. The problem of thermal stresses in irregular cylinders 
is formulated in a manner that permits solution by the use of finite- 
difference methods. The contour integrals of reference 5 expressing 
the single -valued character of the displacements for arbitrary circuits 
around the internal boundaries are written in forms suitable for 
numerical methods of differentiation and integration. The bo undar y 
conditions based on the assumption of force-free boundaries and single- 
valued displacements are formulated in terms of derivatives of the 
stress function as suggested for uniform-temperature problems in 
reference 5, and stress functions are set up in a manner that is an 
extension to the thermal-stress problem of the work of reference 6 
on doubly connected- domains at uniform temperatures. The relaxation 
techniques of reference 7 are used in solving the finite -difference 
problem of determining the stress functions. Details of the method 
are illustrated by examples . The method is applied to a symmetrically 
heated, hollow circular cylinder and also to a hollow circular cylinder 
■with asymmetrical heating to show that the relaxation technique gives 
approximations to exact answers obtaired by direct mathematical methods. 
Comparison of stresses calculated by .he relaxation tec hni ques with 
those determined by exact methods is also made to compare the r elat ive 
accuracy of several numerical method.j of differentiating the stress 
function to obtain stresses . The calculations for the concentric 
cylinder are described in detail sufficient to permit the method to be 
applied to more irregular cylinders. 


SYMBOLS 

The following symbols are used in this report: 
a il' a i2 ;a i3 constants of integration 

modulus of elasticity in tension and compression 

direction cosines of normal drawn outward from region 
bour.ted by plane curve 
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w 


n 

Q 


r >6 


s 

T 

u ,v,¥ 
x,y,z 
a 

^"xy^*xz J ^yz 

Sc> € y ,C z 


a ,0 .a 
x* y* z 


xy 


0 

03 


distance in xy-plane normal to plane curve 
residual in relaxation calculation 
polar coordinates 

arc length of plane curve in xy-plane 

temperature above initial stress -free state 

components of displacement 

rectangular coordinates 

coefficient of linear thermal expansion 

shearing strain components in rectangular coordinates 

unit elongations (strains) in x- } j-, and z -directions, 
respectively 

Poisson's ratio 

normal components of stress parallel to x~, y-, and 
z-axes 

shearing-stress component in rectangular coordinates 

Airy's stress function 

component of rotation about z-axis 


harmonic operator 


Sx 2 Sy 2 


biharmonic operator 


dx* 


-v 2 2 

ox y 


ay 4 


THEORY 

Assumptions 

The following conditions are assumed to exist: 

1. Steady-state heat flow exists with temperatures uniform along 
any line initially parallel to the axis of the cylinder. 



4 


NACA W 2434 


2. Heat sources and sinks are assumed to "be distributed on the 
external and internal boundaries . 

3. The temperature distribution is assumed to be determined by the 
boundary temperatures and Laplace's equation. 

4. The material behaves in an elastic manner. 

5. The variation of the elastic constants (modulus of elasticity., 
Poisson's ratio , and coefficient of thermal expansion) with temperature 
may be neglected in determining the thermal stresses. 

6. Plane sections initially normal to the axis of the cylinder 
remain plane. 

7. The strains and rotation are constant along any line initially 
parallel to the axis of the cylinder. 

The extent to which a particular structure would fulfill these 
conditions would depend on the particular circumstances. The last two 
assumptions are appropriate when the cylinder is long in comparison with 
its cross-sectional dimensions or when end conditions impose suitable 
restraint. The sixth assumption allows bending of the cylinder about 
axes perpendicular to the axial direction of the cylinder in a manner 
that might vary in the axial direction. The last assumption permits 
planes initially perpendicular to the axis of the cylinder to take on 
a warped shape as a result of the deformation but restricts the bending 
to a circular arcj that is, the radius of curvature of lines initially 
parallel to the axis of the cylinder does not vary along the length of 
these lines. The assumption of no variation of rotation in the axial 
direction is equivalent to assuming that the cylinder does not twist. 


Basic Equations 

The case of plane strain where body forces are given by the 
gradient of a potential is treated in reference 5. As is shown by the 
detailed derivation in appendix A, the governing equation of refer- 
ence 5 applies to the thermal-stress problem defined by the preceding 
assumptions. This equation is 


V 4 0 => 0 

where ft is Airy's stress function defined by 


( 1 ) 
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0 x ~ 2 

ay 


a = ^ 
7 Sx 2 


xy 


dxSy 


( 2 ) 


As shown in appendix B, the axial stresses can he calculated from 
the equation 


a z = E(ax + by + c) + V (a x + a ) - a ET 


( 3 ) 


Boundary Conditions 

Two basically different types of auxiliary condition are established 
for the set of physical conditions assumed to exist at the bounding 
cylindrical surfaces. These conditions are: 

(1) Conditions stemming from the nature of the forces applied to 

the surfaces 

(2) Conditions stipulating that the displacements be single-valued 

In calculating the stresses due to the temperature distribution, 
other stresses such as those due to centrifugal force and fluid pres- 
sure are to be calculated separately and the total stresses are to be 
obtained by superposition. All internal and external boundaries are 
therefore postulated to be free from applied forces, and the boundary 
conditions (reference 8, p. 21) become 

a x 1 + T xy m 

T xy 1 + a y 111 



where 


Z 


= 

ds 


dx 

ds 


m 
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Substitution from equations (2) into equations (4) yields 

d 2 {^ dy + 52 ^ dx 

Sy 2 ds 5x3y ds 

£x = 0 

SxSy ds dx 2 ds 

from which 

d d0 
ds 5y 

_d_ ££ 
ds 3x 

along the boundaries . 



Conditions for Single -Valued Displacements 
The defining equations for strains and rotation are 



(5) 


From equations (6), 
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du 

dx 


Sc 


dv 

dx 


2 r *y 


+ 05 


Su = l r 

dy 2 r xy 


05 


dv 

oy 


e 


y 






cto _ 1 °^xy ^ € x 
dx 2 dx dy 

oo5 1 foxy + ^y 

dy 2 dy dx ' 


(7) 


The change in rotation for an arbitrary path P i of integration 
starting at some point (xqjYq) and returning to the same point after 
enclosing (and only enclosing) the internal boundary Ch is 




Oy 


( 8 ) 


is 


From equations (7) and (8), the condition that 05 be single valued 


05 


f.l ^7~xy ^Sc 
2 ox dy. 


i dy- Be, 
dx + U i —^L ^ —3 


2 dy dx/ 


Oyf = o (9) 


By use of equations (A8) and (A9) of appendix A, equation (9) may 
be written 


£ 


ds = - 
<J On 1 




dr 

dn 


ds 


"(10) 
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The change in the x~coraponent of displacement for an arbitrary path 
of integration starting at some point (xq, y Q ) and returning to the 
same point after enclosing the internal boundary is 


is 



( 11 ) 


From equations (7) and (ll), the condition for single-valued u 



«x 


dx + 


xy 


~) 


dy = 0 


( 12 ) 


For the term involving rotation in equation (12), integration by 
parts gives 



■where 

Ho = 0 

because of the single-valued character of u> achieved by imposing 
equation (10) . Furthermore, 
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Equation (12) can now be written 

n r 

% 


Q> < 


€ x<ax + -rxy<3y + -y 


Sr 


xy 


2 • ^ " 2 ” dx 


dx - 


Sc, 


dy\ 


•X 1 .. ^ ~ ^ 


5«, 


y ^ ta -2 y 8y 


<ay + y -=^ &y 


(13) 


In equation (13 ) , integration of the first two terms by parts gives 

n r 


u 

p. 




(• 


c x dx + 


2 'xy 


W 1° lT 1° 

= x€ x +2 ^xy “ 

1 L Jo L Jo 


f 




Sc, 


x\-r— dx + -r- ay 
vox oy 


+ ^7 dJ ) 


where 


HM yr *i = 0 


because the strains are single valued. Equation (13) becomes 




Sc x Se x 

y *sF“ x -^ 


■)dx (- 


y ^ +y ^. x ^p) 


Sy 


8* ’ * li?/ ^ ° 


(14) 


Substituting from equations (A8) and (A9) of appendix A and 
stipulating that the stresses are single valued give equation (14) as 




ds 


( 15 ) 
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If the integrations of equation (15) are performed along the 
internal boundary C^, then from equations ( 5 ), 



Because the integrand of the preceding integral is an exact differential, 
the integral has the same value (zero) for all reducible paths; hence 
for the arbitrary circuit equation (15) becomes 



Similar reasoning from the single -valued character of the y-component 
of displacement leads to the equation 



Equations (10) , (16), and (17) are similar in form to some equations 
derived by Mindlln in reference 9. The equations of that reference 
apply to Integration paths taken along the boundaries . 


Boundary Constants 

The initial objective is to determine the stress function 0 so as 
to solve the boundary-value problem associated with equation (l) . Numer- 
ical differentiation could then be performed to determine the stress 
components according to equations ( 2) • 


The assumption of stress-free surfaces restilted in equations (5). 
With a^ and a^g as constants of integration, these boundary con- 
ditions may be integrated to give along the boundaries 




= a. 


11 


ae5 

a?' ai2 



J 


( 18 ) 
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In general, 

<10 = 2 ^. dx + dy 
dx dy 

and -where is a constant of integration, integration of the pre- 

ceding equation along a boundary and use of equations (18) give 

0i = a^x + a i2 y + a^ (19) 


along a boundary . 

Differential equation (l) and boundary conditions, expressible by 
equations (18) and (19), do not completely determine the stress distri- 
bution. No temperature terms are present, and the values of a^--, , 
a 12 , and a^ are still to be specified. The initial step in obtain- 
ing the compatibility condition (equation (l)) was to raise the order 
of the first differential equation of equations (A2) by differentiating 
once with respect to x and once with- respect to y. These differ- 
entiations, although useful in simplifying the form of the equation, 
require that additional factors be considered. In the case of 
multiply connected regions, the restating equation (equation ( 1 )) does 
not preclude the occurrence of stress distributions due to mechanical 
dislocations (reference 9). The occurrence of mechanical dislocations 
is eliminated by imposing conditions that rotation and displacements 

be single valued. These conditions are to be used In evaluating the 

constants of equation (19). The physical condition's of the problem will 
then be satisfied, for the stress function has been defined so as to 
satisfy the equilibrium conditions j the compatibility condition is 
satisfied by the use of biharmonic functions j the assumption of force- 
free boundaries is satisfied by conditions (18) and (19) on the stress 
function? and the constants a^ , a^ 2 , and a .^3 are to be evaluated 
so as to satisfy the conditions of single-valued rotation and 
displacements. That adjustment of the values of a ^ , and a i3 

is sufficient to satisfy the conditions of single-valued rotation and 
displacements will become evident In the next section. A method of 
determining the constants that will be appropriate for numerical 
techniques is needed. 


Determination of Boundary Constants by Formation of 
Special Solutions 

A method of using special solutions to determine a ^ , a i2 , and 
a ^3 for a doubly connected domain at uniform temperature was described 
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by Prager (reference 6 } . A similar method was suggested for the multiply 
connected domain by Southwell (reference 10). The method of Prager is 
here extended to domains with more than one hole and with temperature 
distributions present. The method is then modified to a form that is 
suitable for numerical techniques . 


Let = 1, 2, . . . k, . . . n; j = 1, 2, 3) be special 

solutions of equation (l) that are defined over the domain bounded by 
the external contour Cq and the internal contours C-]_> Cg* • • • 




■'n* 


Because equation (l) is linear,, the products of the 


and the arbitrary constants a^j may be superposed to give the 
complete solution of the boundary- value problem according to the scheme 


-tt 

i -1 j=l 


^LJ 


( 20 ) 


provided that the boundary conditions for the 0 jM are properly 
selected and that the values of the boundary constants a-jj are prop- 
erly evaluated. 

The boundary conditions for the 0£j must be selected so that the 
function given by equation ( 20 ) will satisfy equations (18) on all the 
boundaries. Because the stresses are given by second derivatives of (, 
the addition of a linear function of the coordinates to 0 will leave 
the stresses unaltered. The assumption is now made that this addition 
is accomplished so that the boundary constants of equations (18) and 
(19) are zero on the external boundary. The boundary conditions on 
the j for the external contour Cq are therefore taken as 

0ij “ - 5 ^ = 0 on C 0 


The boundary conditions for the 03, j on the internal boundaries must 
now be selected so that 3n linearly independent solutions for the 
03 , j will be obtained and so that the function defined by equation ( 20 ) 
will satisfy equations (18) on all the internal boundaries. This 
selection is accomplished as follows: 


¥lj 

Sn 


Let 


hi 


0 


( 21 ) 
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on a] 1 boundaries except 

C k* 

On C k , let 

- 

1! 

X 

S 0kl = dy 
dn ds 

(22) 

0kZ = 

y 

S 0k2 dx 

dn - “ ds 

(23) 

ii 

l 

^-0 

(24) 


Certain geometrical aspects of the choices expressed "by equa- 
tions (21) to (24) are now mentioned. The observation is made, for 
example, that the equation 

dffkl = dy 
dn ds 

is a restriction on in addition to the restriction of = x 

on C k , inasmuch as 0^-1 = x prescribes values of ^ only on the 
line defining the boundary C^, whereas the normal derivative 

dftkl _ dy 
dn ds 

specifies the rate of change of as the boundary is crossed in a 

direction normal to the boundary. 

The slope of the plane 0^ = x in a direction normal to the 
contour is d$yi /on; but because for this plane 0 Tri = x, the 

slope may be written as dx/dn. In the direction n, the slope of 
the surface is given by equation (22) as 

5 $kl _ dy 
dn , ds 

but in the xy -plane, 

dy = dx 
ds dn 

Therefore, in the direction of n, the slope of the surface (in 

equation (22)) has been taken equal to the slope of the plane = x 

a pri because the intersection C’jj. of the cylinder through C^. with the 
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surface 0£j_ lies in the plane 0^j_ = x, the surface 0y| is tangent- 
to the plane 0 yi = x along C’^. The particular solution a^ 0}^ 
therefore defines a surface tangent to the plane 0 y] =. a >-| x. Equa- 
tions (23) require that 0^2 tangent to the plane 0^2 = ^ along 
the intersection C"^ of the cylinder through with the surface 

0k2* and equations (24) require that 0 ^ 3 be tangent to the plane 
0 kg = 1 along the intersection C"‘k of the cylinder through 
with the surface 0 ^ 3 . 

The three special solutions, a kl 0ki, a^J^* ^ a k3^k3> ass °“ 

elated with C k are seen to be tangent to three planes, one of which 

passe's through the y-axis, one of which passes through the x-axis, 
and one of which is parallel to the xy-plane. Determination of a kl , 

a k 2 > and a^g is seen to be equivalent to dete rmining the slope of 

the plane through the y-axis, the slope of the plane through the 
x-axis, and the height of the plane parallel to the xy-plane. Super- 
position according to 

^k = a kl^kl + ®k2^k2 + ^3^3 

is thus seen to satisfy the requirements of equations (18) and (19), 
but the constants a k1 , a k2 > and a^g must be properly chosen. Super- 
position of all the a 1 j0 i j will still leave equations (18) and (19) 
satisfied on C k because of the requirements laid down on all the 0-jj 
by equation (21). 

The 0^ were defined as special solutions of equation (l) , which 
is equivalent to writing 

V 4 0ij => 0 

(i = 1, 2, . . . k, . . . nj j = 1, 2, 3) (25) 

Equations (25) together with boundary-condition equations (21) to 
(24) constitute 3n boundary -value problems for the 3n particular 
solutions 0i j . These individual boundary-value problems with 0^ j 
and 50-y/cto. specified on every boundary are now to be solved by methods 
already described by Fox and Southwell (reference 7). With the assumption 
that the 0^^ are so evaluated, the next step is to calculate the values 

of & ± y 

The method for determining the a^j may be symbolically expressed 
by substituting the complete integral as expressed by equation (20) into 
the contour integrals of equations (10), (16), and (17) . This 
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substitution requires that, on each internal boundary, contour 

integrals involving all the 0-j_ j be formed as coefficients of the ajy . 

Formation of these contour integrals then permits writing simultaneous 
equations for each k'kh boundary (k = 1, 2, . . . n) so that 3n 
equations are obtained for the a-y where the contour integrals 
involving the become coefficients of the ay : 



Because there are exactly 3n equations with exactly 3n unknowns 
(the ay ) , the determination of these constants is sufficient to insure 

the occurrence of single-valued rotation and displacements. 


NUMERICAL TECHNIQUE 

The purposes of this section are (l) to illustrate, by a concrete 
example, the detailed steps by which a solution of a problem may be 
obtained, (2) to show, by comparison with a problem for which the 
solution is already known, that the application of relaxation procedures 
to the finite-difference method outlined in the preceding part of this 
report yields a method giving results that approximate the correct 
answers, and (3) to present the results of an investigation of some 
factors affecting the accuracy of the answers for some particular 
conditions encountered In the relaxation solution of a thermal stress 
problem. 
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Illustrative Example 

The detailed instructions enabling the relaxation calculation of 
thermal stresses in an irregular cylinder are presented in appendix C. 
The particular problem illustrated is that of a concentric circular 
cylinder possessing an asymmetrical temperature distribution; however, 
the description of method may be applied to cylinders of more complex 
shape. The choice of the concentric circular cylinder enabled the 
relaxation work for the stress functions to be confined to a 90° sector. 
There would be no fundamental distinction in carrying out the cal- 
culations for an irregular profile - relaxation would merely have to be 
performed over the entire cross section entailing more labor. (Although 
circular boundaries were involved, the advantages of using polar 
coordinates were not utilized. The use of rectangular coordinates in 
the presence of circular boundaries involves boundary technique problems 
typical of a more irregular region. ) 

The results of the relaxation calculation (reference 11) according 
to Laplace's equation for the temperature distribution are presented 
in figure 1. Results of relaxation calculations for special solutions 
of the biharmonic equation are presented in figures 2 to 4. Contour 
integrals were calculated as illustrated by table I and the stress 
function is presented In table H. Tangential stresses (table II) 
were calculated by computing the second derivatives of the stresB func- 
tion with respect to radius (reference 8, p. 53) according to the 
5-point method of reference 12 described in appendix D. 

Exact values of radial and tangential stress were calculated as 
indicated in appendix E and are also listed in table II. The maximum 
tangential stress Is seen to be much larger than the maximum radial 
stress. The error in the relaxation calculation of the maximum tan- 
gential stress was about 5 percent. Comparison of the values of tan- 
gential stress calculated by the exact method with values calculated 
by the relaxation method is also presented in figure 5. 


Investigation of Factors Affecting Accuracy 

Several factors influencing the accuracy obtainable in calculating 
thermal stresses by the method just described were investigated by 
applying various calculation techniques to a problem for which answers 
could be calculated by exact mathematical methods. The example chosen 
consisted of a concentrio cylinder subjected to a symmetrical tem- 
perature distribution. Details of the relaxation calculation are 
presented in appendix F. 

The results of the relaxation calculation for the temperature dis- 
tribution are presented in figure 6, Because dimensions of the oylinder 
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were exactly the same as those of the illustrative example previously 
discussed, the special solutions $H) $12 > 8®d 013 determined in that 

example are usable in the present example. This situation illustrates 
an Important feature of the use of relaxation methods in calculating 
thermal stresses; that is, once the time consuming biharmonic relaxation 
work for the special functions $±i, $±2t 0i3 ^ ias ^ een completed for a 
given shape of body, relatively little extra work is required to study 
the effects on ther mal stresses of ohanges in temperature distribution. 

For the symmetrical temperature distribution of figure 6, the 
boundary constants a-| -\ and a]_2 were found to vanish. The values of 
a 13 for the paths a, b, and o of figure 4 are presented in table III. 
The Airy stress function and the stresses calculated from it by numerical 
differentiation according to appendix D are listed in table IV. 

Exact values of Airy's stress function were calculated according to 
appendix G-. The value of the arbitrary constant D in equation (G4) 
was adjusted to give 0 = 0 for r = 12. Results are listed in 
table IV. Second derivatives of the exact values of Airy's BtresB func- 
tion were then calculated by the numerical methods of appendix D to give 
the tangential stress values listed in table IF. Exact values of tan- 
gential stress were calculated using the second of equations (El) and 
are also listed in table IF. 

Comparison of the errors in tangential stress of table IF shows that 
the errors associated with the numerical differentiation of the relaxa- 
tion calculated Airy function were much larger than those associated with 
the numerical differentiation of the exact Airy function. (Errors In the 
mrimuTn stress were 21.5 and 6.5 percent, respectively.) This comparison 
suggests that the relaxation calculated Airy function was an important 
source of error. 


The relaxation calculated stress function was calculated as the 
product aj _3 023 (appendix F) . An exact value of aj _3 was calculated 
by observing that the exact values of 0 in table IF range from zero 
at r = 12 to -467,627 at r = 4; whereas in figure 4 the values of 
013 range from zero at r = 12 to 1000 at r = 4. The exact value 
of a 23 is therefore 


a 13 


-467.627 

1000 


-467.627 


As Ind icated by the data of table III, the errors associated with 
the processes of integrating difference quotients to calculate the 
boundary constants can be significant but are small if averaged over 
several paths. 
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The extent to -which errors in the relaxation calculation of the 
Airy function influence errors in the stresses was evaluated hy ro und - 
ing off the previously determined exact values of Airy's stress func- 
tion to three significant figures. Corresponding stresses and errors 
are presented in table 17. In general, the errors associated with the 
three-figure Airy function are seen to he significantly lower than those 
associated with the relaxation calculated Airy function. Apparently, 
the equivalent of three -figure accuraoy was not achieved in the relaxa- 
tion calculation of the Airy function. Comparison of the values of 
tangential stress calculated hy the exact method with values calculated 
hy the relaxation method is also presented in figure 7 . Improved 
accuracy could he accomplished hy (1) further reduction of residuals 
with the introduction of another significant figure, (2) the use of 
a finer net spacing, (3) the use of more elegant finite -difference 
methods, or (4) scane combination of (1), (2), and (3). A critical 
discussion of some factors influencing the aocuracy of relaxation pro- 
cedures is contained in reference 13. 


CONCLUDING REMARKS 

This investigation has yielded a numerical method for calculating 
thermal stress In a cooled irregular cylinder possessing one or more 
cooling passages under steady-state temperature conditions. Application 
of the method to structures such as internally cooled turbine blades is 
suggested. The use of relaxation methods and elementary methods of 
finite differences has been found to give approximations to oorrect values 
when oompared with previously known solutions for concentric circular 
cylinders possessing symmetrical and asymmetrical temperature distri- 
butions. 


Lewis Flight Propulsion Laboratory, 

National Advisory Committee for Aeronautics, ' 
Cleveland, Ohio, May 10, 1951. 
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APPENDIX A 

DERIVATION OF BIHARMONIC EQUATION 

An outline of the derivation of the governing partial differential 
equation for the stress function is given in reference 4. A detailed 
derivation using the conditions imposed in the section on assumptions 
follows : 


are 


The defining equations for the normal strains (reference 8, p. 7) 


£ _ Bu 

X "S 




€y By 


$ > 


By 

By 

Bw 


£z Bz 


(Al) 


and the defining equations for the shearing strains are 

7xy By Bx 


_ Bu Bw 
7xz ~ B^ + Bi 


7yz 


By By 
Bz By 


(A2) 




The defining equation for the rotation (reference 8, p. 162) is 


The assumption that the strains and rotation are constant in the 
direction of the z-axis (axis of the cylinder) permits writing 



20 


NACA TN 2434 


The conditions that the strains "be compatible with displacements 
specified by u, v, and w are (reference 8, p. 196) 


By 2 Bx 2 dxd y 

+ fi- fZz* 

Bz 2 By 2 ByBz 

B^z + jffx = a2 7x Z 
5x 2 Bz 2 BxBz 

^fx B l±f. ^zz + + 5y xyj 

ByBz 2 Bx \ Bx By Bz / 

B 2g y = 1 B_ /^yz _ _Jrxz + W ) 
BxBz 2 By \ Bx By Bz / 

a2g z = 1 B /By yz + f£xz _ ^x^j 
BxBy 2 Bz \ Bx By Bz / 




> 


J 


From equations (A3) and (A4) , equations (A5) beoome 

a^x b 2 ^ y b 2 7 ^ 

By 2 Bx 2 Bxdy 



B 2 e 


Bx 


o 


B Bw 
3x 


5 0 


0 


B B^> 
By Bz 


o 


B g e z 

BxBy 


0 


> 


(A5) 


(AS) 
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The first of equations (A6) remains whereas the fourth and fifth of 
equations (A6) vanish identically. The second, third, and sixth of 
equations (A6) may be written 

3^w _ d^w _ S^w n /» 7 \ 

Sy 2 Sz ” 3x 2 dz " = 

Because of the assumption that plana sections initially normal to 
the axis of the cylinder remain plane after the deformation, the 
displacement in the z-directicn may be Written 


w = f (z) [jax + by + c] 


showing that equations (A7) vanish identically. The preceding proof 
that the last five of equations (A6) vanish identically is essentially 
a demonstration that these compatibility conditions are satisfied as 
a consequence of the initial assumptions. Of equations (A6), the 
only nonvanishing equation is now the compatibility condition 


d 2e x S2e y & 7 X y 


(A8) 


The generalized Hooke’s law equations (reference 8, p. 204) are 


*x 

■i[ 

°x " u (°y 
> 

+ a z )] + aT 

€ , r 

-if 

a - u (a 

+ 0)1 + aT 

y 

E L 

y ' x 

z J 

e 

-if 

o - u (a 

+ 0 )"| + aT 

z 

E l 

z ' X 

y J 


J 


(A9) 


Elimination of O z between these equations yields 
e x = i [(1 - 1?2) o x - h (1 + V) a y + (1 + V) oET - UE cj 

Cy * i [(1 - v2) O y - v (1 + v) O x + (1 + V) OET - WE c z ] 


(A10) 


Airy's stress function is defined by 
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3y2 

*>■3 




> 


_ - 

BxBj 


( 2 ) 


Equations (AlO) and (2) are used with the direct method of calculat- 
ing thermal stresses; that is, the equilibrium, boundary, and compati- 
bility conditions are used without regarding temperature terms aB body 
and surface forces. The equilibrium equations (reference 8, p. 195) are 
therefore written without body forces . These equations were reduced to 
those of reference 8, page 21, by using the assumption that plane sections 
initially normal to the axis of the cylinder remain plane after the 
deformation, so that 


d°x . d T xy 

-& +_ 3T ■ 


0 


= 0 

By Bx 


(All) 


Substitution of expressions (2) in equations (All) shows that the 
stress function has been defined so as to satisfy the equilibrium equa- 
tions identically. The condition remaining to be satisfied by is 
the compatibility condition (A8). Substitution of expressions (2) in 
equations (AlO) yields 

“0 - ^ - UE£ *) + a + 

‘y » | - »|^| - + (1 + V>) OT 

The shearing strain is expressed (reference 8, p. lo) in terms of 
the shearing stress by 



* 


_ ,2(1 f?± 

7 xy E xy 


2(l+u) B2jzi 
E BxBy 


(A13) 


CRT? 



2193 


MACA OT 2434 


23 


By use of equations (A6), (A12), and (A13), equation (A8) can be 
written 

V 4 # = - _SJL 7 2 T 

1 - u 

For the assumed tempera 1 ture conditions, 

V 2 T = 0 

and hence the governing equation for $ is 

V 4 0 = 0 (1) 
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APPEKDIZ B 


CALCULATION OF AXIAL SERESS 


After O x and Oy have been calculated according to equation (2), 
the normal stress O z in the axial direction may be determined. The 
third of equations (A9) is 

e z = I [a z - U (a* + o 7 )] + csr 

from, which 

a z = € z E + v (o x + O y ) - a ET (Bl) 


As stated in the section Assumptions , the cylinder is free from 
applied forces and therefore the total force in the axial direction 
and the x- and y- components of a bending couple are zero. 



dx dy = 0 


dx dy e 0 


dx dy = 0 


> 


(B2) 


Because of the assumption that plane sections initially normal to 
the axis of the cylinder remain plane after the deformation, the 
strain in the z-direction may be written 

e z = ax + by + c (B3) 

where a, b, and c are constants. 

Substitution of expressions given by equations (Bl) and (B3) in 
equations (B2’) gives 
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a jT Itady + l ’ fl y dx dy + c ^ J dx dy 


| »( a z + °j) - a Bt] dx dy 


7/ 


a II x 2 dx dy + 1) 

'A ^°A 


// Iy4iay + °// 


c / / x dx dy 
A 


= - | /jf x J»(a x +Oy) - ct EtJ dx 


4y 


r (B4) 


. Jx y dx dy + Id y 2 dx dy + c 


y dx dy 


[“ (°x +°y) - 


dx dy 


Equations (B4) are tliree simultaneous equations from which. the 
constants a^ b^ and c may he evaluated. Substitution of the 
expression for € z given in equation (B3) for c ^ ^- n - QOus-'ki 01 ! (Bl) 

gives 

a z = E(ax + by + c) + V (a z +a y ) - a EE (3) 

from which the axial stresses may be determined. 
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APPENDIX C 

DETAILS OF HELAXATION CALCULATION OF TANGENTIAL STEESS 
IN CONCENTRIC CYLINDER WITH ASYMMETRIC 
TEMPERATURE DISTRIBUTION 

The temperature distribution -was assumed to be determined by tem- 
peratures in °F of zero on the internal boundary of 4- inch radius and 
500 + 1000 cos 0 on the external boundary of 12- inch radius. The 
values cf the elastic constants were assumed to be 

a = 8.0 X 10-6 (in./in. /°F) 

E = 17.5 X 10 s (lb/sq in.) 

V = 0.3 

The results of the relaxation calculation for the temperature 
distribution according to the technique described in reference 11 is 
presented in figure 1. In all calculations, the origin of coordinates 
was located at the center of the cylinder. In the present case (one 
internal boundary), the indices of equations (21), (22), (23), and 
(24) became 


i = k = n = 1 
J = 1, 2, 3 

In order to eliminate the use of decimals in the work of relaxation, 
equations (21), (22), (23), and (24) were modified as follows: 

On the exterior boundary, 


on 


on 


On the interior boundary, 

011 = 250 x 

012 = ^50 y 

0X3 = 1000 


. 350 a* 
on ds 

. -250 52 
on ds 


5 013 

3n 


= 0 
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The boundary value problems just defined for the biharmonic functions 
^11* $12 > an ^ $13 "were then solved by the techniques of references 7 and 
11 to yield the solutions presented in figures 2 to 4. (The dimensions of 
the cylinder were such that the distance between nodal points could be 
conveniently taken as 1 inch. For more general cases a method of handling 
dimensions , nodal distances, and derivatives, such as described in refer- 
ence 7, can be followed.) 

The next step is to calculate the values of the constants a ^ , 
a -| ? , and a -^ according to eqLuations (26), (27), and (28). To that end, 

the appropriate contour integrals were calculated as illustrated for 
functions involving and T along path c in table X. The 

numbering system for points in the tables corresponds to the numbering 
system exhibited by the small diagram in figure 2 , where j is the station 
along the path of integration and j = 1 on the positive x-axis. 

With the use of averages of contour integrals for paths a, b, 
and c, the simultaneous equations (26), (27), and (28) become 

1241.3 a-j^ = -200 X 2853 

8466.3 ai2 «* 0 

8466.3 a-jj_ = -200 X 9532.3 

from which 

a^ = -459.7 

a 12 = 0 
a-j = -225.2 

According to equation (20), 

$ = a ll $11 + a 12 $12 + a 13 $13 

from which 

0 = -225.2 fS -459.7 0 
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The solution of this equation for points along the positive 
x-axis (6 = 0°) is presented in table II. Tangential stresses were 
calculated from the relaxation solution for the stress function by 
taking second derivatives with respect to radius according to the 
formulas of appendix D. These results are also presented in table II. 
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APPENDIX D 

NUMERICAL DIFFERENTIATION OF AIRY’S STRESS FUNCTION 

The purpose of this section is to present the methods used in 
approximating to the second derivatives of Airy's stress function. The 
function to he differentiated is expressed by 

y = f(x) 

■where In the present application y is Airy's stress function and x 
is radial distance on a cross section of the hollow concentric cylinder. 
Where p is an integer, h is a uniform tabular interval of the 
independent variable, and Xq is an arbitrary point from which the 
distance ph is measured, let 

Xp = Xq + ph (Dl) 

The values of second derivatives were calculated by the five- 
point formulas of reference 12. For the various values of p, these 
formulas are; With p = 0, 

D 2 y 0 = _L_ (35 y 0 - 104 y x + 114 y 2 - 56 y 3 + 11 y 4 ) (D2) 

12h 2 

■with p = 1, 

D 2 y! = (11 y 0 - 20 y]_ + 6 y 2 + 4 y 3 - y 4 ) (D3) 

12b~ 

with p = 2, 

D 2 y 2 = (- y 0 + 16 y x - 30 y 2 + 16 y 3 - y 4 ) (D4) 

12 n 

with p = 3, 

D 2 y 3 = ■ (- y Q + 4 y 2 + 6 y 2 - 20 y 3 + 11 y 4 ) (D5) 

12 h. 

vlth P = 4, 

D 2 ^4 = ~~2 Cl1 7 0 - 56 y x + 114 y 2 - 104 y 3 + 35 y 4 ) (D6) 

12h 

Equation (D2) was used to calculate derivatives at r = 4, equa- 
tion (D3) was used for r = 5, equation (D4) was used for r = 6 to 10, 
equation (D5) was used for r = 11, and equation (D6) was used 
for r = 12. 
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APPENDIX E 


EXACT DETERMINA.T ION OF STRESSES FOR CONCENTRIC CYLINDER 

WITH ASYMMETRIC TEMPERATURE DISTRIBUTION 

The radial and tangential stresses , respectively, in the case of a 
symmetrical temperature distribution are given in reference 8 on 
page 372 by 


°r - 


CtETa 


- 


2(1 - U)ln - 

£L 


clETj 


- In - - 
r 


yZ a 2 
b - a 


( ' rj ln » 


£ 


2(1 -U)ln| 


ln — - 


2 

b - a 




1 

( 1 + ^) ln ! 


r w 


■where a is the radius of the internal boundary and b is the radius 
of the external boundary, and where T^ is the excess of the inner 
boundary temperature over the temperature of the external boundary. 

From reference 4, the stresses in a concentric circular cylinder 
possessing the temperature distribution 


T = A ln r 


w 

+ A o + zI 

n=l 


(A n r 31 + C n r“ n ) cos nS + 


(- B n r 11 + D n r“ n ) sin n8 
where A, B, C, and D are constants, are 


(E2) 
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b 2 - r 2 
b 2 - a 2 


fr 2 - a 2 ' 


(a 2 + b 2 )r3 


sin 0) 


b 2 + r 2 
b 2 - a 2 


f4a 2 r 2 -(5r 2 - b 2 )(r 2 + a 2 )l ( c cos Q + D 

»> L (a 2 + b 2 )^ J 



The temperature distribution in °F is specified by T = 0 on the 
internal boundary and by T = 500 + 1000 cos 0 on the external boundary. 
The stresses -will be calculated by superposing stresses for a symmetrical 
temperature distribution on those calculated for a special asymmetrical 
distribution. The given temperature distribution is resolved into two 
components: a symmetrical component defined by T a 0 on the interior 

boundary and T = 500 on the exterior boundary? and the special asym- 
metrical component defined by boundary temperatures of T = 0 on the 
interior boundary and T =* 100C cos 9 on the exterior boundary. 

The stresses calculated for the symmetrical component of the 
temperature distribution were calculated according to equations (El) . 

Now equation (E2) is to be written so as to satisfy tbe boundary 
co nd itions on the asymmetrical component of the temperature distri bution. 
The boundary conditions are 

at r = a T = 0 

at r = b T = 1000 cos 9 

In carder to satisfy these conditions and equation (E2), let 

n = 1 

A = Aq = B n = D n = 0 


(£4) 
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Equation (E2) is thereby reduced to 


On r = a, 



from which 


Ai = 


£1 

a 2 


On r = b, 


from which 


■ 1000 cos 8 = 



A}*) 



1000 


Hi b + Hi = 1000 

a 2 b 


or 

C, b2 .) = 1000 

a 2 b 

from which 

q 1000 a 2 b 

1 b 2 - a 2 


Use of equations (E4) and of the preceding value of C-j_ 
tions (E3) shows that the components of stress due only to the 
perature distribution specified by T = 1000 ccs 9 on r = b 


in equa- 
tem- 
are 
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= 1000 a. 2 !) (b 2 - r 2 )(r 2 - a 2 ) cos Q 

r 2(1 - b ) r3 b 4 - a^ 


1000 aE a 2 b 4a,2r2 - (3r^ - b^)(r2 + a^) 
° d ~ 2(1 -b) r3 b 4 - a 4 


► 

QOS 9 

J 


(E5) 


Radial and tangential stresses vere superposed from stresses cal- 
culated according to equations (El) and (E5) to give the exact stresses 
listed in table II. 
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APPENDIX F 

DETAILS OF RELAXATION CALCULATION OF TANGENTIAL STRESS 
IN CONCENTRIC CYLINDER WITH SYMMETRICAL 
TEMPERATURE DISTRIBUTION 

The temperature distribution was assumed to he determined by tem- 
peratures of 0° F on the internal boundary of 4- inch radius and 500° F 
on the external boundary of 12-inch radius. The relaxation solution for 
the temperature distribution is given in figure 6. The values of the 
elastic constants were assumed to be the same as those used in 
appendix C, namely , 


a = 8.0 X 10“ 6 (in./ in./°F) 

E = 17.5 x 10 6 (lb/sq in.) 
v = 0.3 

As in the asymmetrical case treated in appendix C , the indices of 
equations (21) to (24) were taken as 


i = k = n = 1 
J « 1, 2/3 

and these equations were modified as follows: 
On the exterior boundary, 


011 = ^12 = ^13 
On the interior boundary, 

011 = 250 x 
/■}_2 = 250 y 
013 = iOOC 


S/ll S 012 S 013 


Sn 


Sn 


S/n 

Sn 

S 012 

dii 

^013 

Sn 


Sn 


350 3Z 
ds 


>250 ~ 
ds 


= 0 
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These boundary values are the same as those used for the preceding 
problem involving the asymmetrical temperature distribution. The 
biharmonic functions $ 12 * an & 013 are Presented in figures 2 

to 4. Calculation of the contour integrals by the methods described in 
appendix C and solution of equations (26) to (28) show that 

a ll = a 12 = 0 

and that the values of a^ for paths a ; b, and c are as given in 
table III. The values of Airy's stress function as determined' entirely 
by the numerical method were then calculated using the relaxation values 
of 0^3 and the average of for paths a, b 3 and c according to 

0 = a 13 013 


Values so calculated are presented in table IV. Tangential stresses 
were calculated from the relaxation solution for the stress function 
by taking second derivatives with respect to radius according to the 
formulas of appendix D. These results are also presented in table IV. 
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APPENDIX G 


EXACT DETERMINATION OF AIRY’S STRESS FUNCTION FOR 
CONCENTRIC CYLINDER WITH SYMMETRICAL 
TEMPERATURE DISTRIBUTION 


The purpose of this section is to present the formula used in 
calculating Airy 's ' stress function for the concentric circular cylinder 
with a symmetrical temperature distribution. A possible form for the 
stress function is given in reference 8 on page 55 as 

f) = A In r + Br 2 In r + Cr 2 + D (Gl) 

where A, B, C, and D are constants and the corresponding stress 
components are given by 

O r = — + B(l + 2 In r) + 2C 
r 2 

a e = - — + B(3 + 2 In r) + 2C 
r 2 



The thermal stresses are given by equations (El) , For comparison 
with the first of equations (G2 ) , the first of equations (El) Is written 

/■ 

aETi a 2 b 2 


a r =<• 


In ^ 


|2(1 - u)ln — b 2 - a 2 a [ r 2 

St 


aETi 


2(1 - v)ln£ 2 


aETj 


2(1 - V)ln | 


i - (1 + 2 In r) + 


i p2 

| ■ In b ■ ^ 5 


ln ! 


(G3) 


Comparison of equation (G3) with the first of equations (G2) shows 
that necessary conditions on the constants of equation (Gl) are 
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aSTj a 2 b 2 

2(1 - U)ln - b2 " a2 

£L 



aET.» 

B = ± 

2(1 - U)ln | 


1 

2 


C = 


aETj 


2(1 - W)lnf 


1 (-i 

2 \ 2 


In Id 



Because the stresses are calculated from, derivatives of Airy’s 
stress function, the value of the constant D in equation (Gl) 
may be taken arbitrarily. The values of A, B, and C as just deter- 
mined are substituted in equation (Gl) to obtain 




oETj 


2(1 - u)m| 


-|r2 inS 


— r 2 


- 


a° | 

V 2 


l 2 

b 2 - a 2 


- b 2 In r 


) 




/ + ID 


(G4) 
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3 ) 

-i 


TABLE I - COHTOUH IHTH3RALS ALOB3 PATH c 


(a) Integrala lnvolrlng 7 2 $;q. 
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v z 6 x 
0 2 


v 2 |5 i 
n+ 5 
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X 

a . v 2 jS , 


y 

dv 2 |jj 


da 

dn 

z de 

y da 

dn 

1 

-113 

-113 

0 

-96 

-170 

74 

5.S 

0 

407.0 

0 

0 

0 

2 

-103 

-113 

10 

-81 

-105 

24 

5.5 

55*0 

132.0 

i 

10.0 

24.0 

3 

-106 

-103 

-3 
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-142 

59 

5.5 
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-127 
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-51 

-92 

41 
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225.5 
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164.0 
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-29 
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-27 

-59 

32 

5 

115.0 

160.0 

- 5.5 

- 126.5 

- 176.0 

40 

-57 

-29 

-28 

-26 

-59 

33 

5.5 

- 154.0 

181.5 

-5 

140.0 

- 165.0 

41 

-85 

-57 

-28 

-51 

-92 

41 

5.5 

- 154.0 

225.5 

-4 

112.0 

-^ 64.0 

42 

-108 

-85 

-21 

-71 

-127 

56 

5.5 

- 115.5 

308.0 

-3 

63.0 

- 169.0 

43 

-103 

-106 

3 

-83 

-142 

59 

5.5 

16.5 

324.5 

-2 

- 6.0 

- 118.0 

44 

-113 

-103 

-10 

-81 

-105 

24 

5.5 

- 55.0 

132.0 

-1 

10.0 

- 24.0 







0 


0 

7540.0 


1914.0 

0 



Total 
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TABL2 I - COIWOUR IOT20RAIS ALOIKJ PATH c - CONCLUDED 


(b) Integrals Involving T. 


4T. ST. 

T i T ! 

as n+g n-^ on 


387 

387 

0 

491 

266 

225 

5.5 

0 

418 

387 

31 

504 

286 

218 

5.5 

170.5 

472 

418 

54 

540 

342 

198 

5.5 

297.0 

539 

472 

67 

S92 

416 

176 

5.5 

368.5 

609 

539 

70 

652 

496 

156 

5.5 

385.0 

676 

609 

67 

714 

573 

141 

5.5 

368.5 

544 

676 

-132 

645 

573 

72 

5 

-660.0 

422 

544 

-122 

522 

438 

84 

4 

-488.0 

315 

422 

-107 

411 

316 

95 

3 

-321.0 

228 

315 

-87 

316 

216 

100 

2 

-174.0 

166 

228 

-62 

239 

143 

96 

1 

-62.0 

123 

166 

-43 

181 

99 

82 

0 

0 

90 

123 

-33 

137 

74 

63 

-1 

33.0 

59 

90 

-31 

98 

53 

45 

-2 

62.0 

22 

59 

-37 

58 

26 

32 

-3 

111.0 

-22 

22 

-44 

13 

-10 

23 

-4 

176.0 

-71 

-22 

-49 

-36 

-54 

18 

-5 

245.0 

-85 

-71 

-14 

-105 

-54 

-51 

-5.5 

77.0 

-94 

-85 


-115 

-67 

-46 

-5.5 

49.5 

-97 

-94 


-121 

-72 

-49 

-5.5 

16.5 

-96 

-97 


-123 

-71 

-52 

-5.5 

-5.5 

-93 

-96 

3 

-122 

-66 

-56 

-5.5 

-16.5 

-93 

-93 

0 

-121 

-63 

-58 

-5.5 

0 

-96 

-93 

-3 

-122 

-66 

-56 

-5.5 

16.5 

-97 

-96 


-123 

-71 

-52 

-5.5 

5.5 

-94 

-97 


-121 

-72 

-49 

-5,5 

-16.5 

-85 

-94 


-115 

-67 

-48 

-5.5 

-49,5 

-71 

-85 

14 

-105 

-54 

-51 

-5.5 

-77.0 

-22 

-71 

49 

-36 

-54 

18 

-5 

-245.0 

22 

-22 

44 

13 

-10 

23 

-4 

-176.0 

59 

22 

37 

58 

26 

32 

-3 

-111.0 

90 

59 

31 

98 

53 

45 

-2 

-62.0 

123 

90 

33 

137 

74 

63 

-1 

-33.0 

166 

123 

43 

181 

99 

82 

0 

0 

228 

166 

62 

239 

143 

96 

1 

62.0 

315 

228 

87 

316 

216 

100 

2 

174.0 

422 

315 

107 

411 

316 

95 

3 

321.0 

544 

422 

122 

522 

438 

84 

4 

488.0 

576 

544 

132 

645 

573 

72 

5 

660.0 

609 

876 

-67 

714 

573 

141 

5.5 

-368.5 

539 

609 

-70 

652 

496 

156 

5.5 

-385.0 

472 

539 

-67 

592 

416 

176 

5.5 

-358. 5 

418 

472 

-54 

540 

342 

198 

5.5 

-297.0 

387 

418 

-31 

504 

286 

218 

5.5 

-170.5 


1237.5 

1199.0 

1089.0 

968.0 

858.0 

775.5 

360.0 

336.0 

285.0 

200.0 

96.0 

0 

-63.0 

-90.0 

-96.0 

-92.0 

-90.0 

280.5 

264.0 

269.5 

286.0 

308.0 

319.0 

308.0 

286.0 

269.5 

264.0 

280.5 
-90.0 
-92.0 
-96.0 
-90.0 
-63.0 

0 

96.0 

200.0 

285.0 

336.0 

360.0 

775.5 

858.0 

968.0 

1089.0 

1199.0 


0 15,843.5 


0 

31.0 

108.0 

201.0 

280.0 

335.0 

-726.0 

-671.0 

-588.5 

-478.5 

-341,0 

-236.5 

-181.5 

-170.5 

-203.5 

-242.0 

-269.5 

-70.0 

-36.0 


-9.0 

-36.0 

-70.0 

-269.5 

-242.0 

-203.5 

-170.5 

-101.5 

-236.5 

-341.0 

-478.5 

-588.5 

-671.0 

-726.0 

335.0 

280.0 

201.0 

108.0 

31.0 


-6527.0 


o o 
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TABLE II - STRESS ALONG POSITIVE X-AXIS IN CONCENTRIC CIRCULAR 
CYLINDER WITH ASYMMETRICAL TEMPERATURE DISTRIBUTION 


Radius 

(in.) 

Stress 

function 

Relaxation 
tangential 
stress 
(lb/sq in.) 

Exact stresses 
(ib/sq in.) 

Error 8 - in 
relaxation 
stress 
(lb/sq in.) 

Radial 

Tangential 

4 

-684,900 

121,012 

0 

126,988 

-5,976 

5 

-652,832 

65,289 

18,126 

69,800 

-4,511 

6 

-554,142 

25,558 

22,172 

35,409 

-9,851 

7 

-428,561 

2,370 

20,982 

11,675 

-9,305 

8 

-299,828 

-11,372 

17,673 

-6,277 

5,095 

9 

-182,297 

-23,815 

13,498 

-20,741 

3,074 

10 

-88,278 

-32,635 

9,013 

-32,929 

-294 

11 

-26,440 

-35,852 

4,473 

-43,542 

-7,690 

12 

0 

-33,617 

0 

-53,012 

-19,395 


Positive signs denote relaxation values larger in absolute 
magnitude than exact values j negative signs denote relaxation 
values smaller in absolute magnitude than exact values. 


TABLE III - VALUES OF a^ AND 
ASSOCIATED ERRORS 



a l3 

Error 

(percent) 

Exact 

-467.627 


Path a 

-555 .0 

18.7 

Path b 

-423.8 

9.4 

Path c 

-416.6 

10.9 

Average of 
a,"b j c 

-465.1 

.5 





CO 


■TABLE IV - NUMERICALLY CALCULATED STRESSES AND ASSOCIATED ERRORS FOR CONCENTRIC 


CYLINDER WITH SYMMETRICAL TEMPERATURE DISTRIBUTION 


Radius 

(in.) 

Exact values 
of Airy's 
stress 
function 

Relaxation 
values of 
Airy's 
stress 
function 

Exact values 
of tangential 
stress 
(lb/sq in.) 

5 -point method 
for tangential 
stress using 
relaxation 
values of 
Airy's stress 
function 
(lb/sq in.) 

5-point method 
for tangential 
stress using 
exact values 
of Airy's 
stress 
function 
(lb/sq in.) 

5 -point method 
for tangential 
stress using 
Airy’s func- 
tion exact to 
three 
figures 
(lb/sq in.) 

Stress 

Error a 

Stress 

Error 8 - 

Stress 

Error 8 

4 

-467,627 

-465,100 

66,988 

81,357 

14,369 

62,620 

-4368 

60,600 

-6,400 

5 

-439,879 

-438,124 

36,582 

42,286 

5,704 

36,989 

407 

36,600 

0 

6 

-374,651 

-367,894 

17,284 

14,844 

-2,440 

17,254 

-30 

17,600 

300 

7 

-291,678 

-281,851 

3,636 

-620 

-4,256 

3,627 

-9 

3,800 

200 

8 

-204,799 

-195,807 

-6,746 

-8,449 

1,703 

-6,766 

20 

-7,300 

600 

9 

-124,507 

-118,135 

-15,058 

-15,969 

911 

-15,054 

-4 

-13,800 

-1,300 

10 

-59,152 

-56,277 

-21,964 

-21,588 

-376 

-21,965 

1 

-23,400 

1,400 

11 

-15,677 

-15,813 

-27,864 

-24,845 

-3,019 

-27,883 

19 

-27,400 

-500 

12 

0 

0 

-33,012 

-25,780 

-7,232 

-32,776 

-236 

-26,400 

-6,600 



Positive quantities denote approximate values larger in absolute ma g nitude than exact val- 
ues j negative quantities denote approximate values smaller in absolute ma gni tude than exact 
values . 
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Figure 1. - Results of relaxation calculation of asjmnetrical temperature distribution. Temperatures in ° 7 : on Inner boundary, 

T » 0; on outer boundary, T = 500 + 1000 cos 0. Temperature -values on integration paths are obtained, by averaging temperature 
values at four nearest nodal points. 
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Figure 2. - Besults of relaxation, calculation of jtf’xi* Values of V jzi 

paths are obtained by averaging V 2 0 values at four nearest nodal points. 

V 2 0 = + ^2 + ^3 + " i f' 
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Figure 3. - Results of relaxation calculation of 0^. Values of V 2 $ on integration 
paths are obtained by averaging V 2 0 values at four nearest nodal points. 

V 2 {? = 01 + lz + 03 + 04 “ 4 0- 
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Radius, In. 

Figure 5. - Comparison of tangential stresses calculated by relaxation method with exact stresses; 
concentric cylinder with asymmetrical temperature distribution. 
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Figure 7. - Comparison of tangential stresses calculated by relaxation method with exact stresses 
concentric cylinder with symmetrical temperature distribution. 



